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We study a mechanism that generates the baryon asymmetry of the Universe dur- 
ing a tachyonic electroweak phase transition. We utihze as sole source of CP violation 
an operator that was recently obtained from the Standard Model by integrating out 
the quarks. 



I. INTRODUCTION 

The requirements for the dynamical creation of the baryon asymmetry in the Universe 
are stated by the Sakharov conditions: Violation of baryon number conservation, violation 
of charge conjugation (C) and charge-parity conjugation (CP) symmetry, and departure 
from equilibrium While there exists a plethora of possible explanations for the baryon 
asymmetry, for a review see e.g. j2|, all of them have to introduce physics beyond the 
Standard Model (SM) in order to provide all these ingredients. Besides, most of them operate 
at energy ranges that are not testable by experiments in the near future, and hence become 
difficult to falsify. In principle, the SM contains all three necessary ingredients mentioned 
before. Baryon number violation is supplied by the weak anomaly |3[ and weak interactions 



violate C maximally and violate CP through the Kobayashi-Maskawa mechanism j^. Finally, 
departure from equilibrium could occur due to expanding Higgs bubbles during a first- 
order electroweak phase transition, and the corresponding mechanism is called electroweak 
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baryogenesis [sl- However, the experimental lower bound on the Higgs mass implies that 
there did not occur a strongly first-order phase transition as required but a crossover in the 
SM ^|. Even if the phase transition would be first-order, the CP violation from the CKM 
matrix is not strong enough in the Higgs bubble walls to generate a substantial baryon 
asymmetry. While this rules out electroweak baryogenesis in the SM, there are extensions 
of the SM which overcome these shortcomings while still remaining close to the SM. This 
includes e.g. the MSSM with CP violation in the chargino [?, 8| or neutralino sector jol or 



singlet extensions of the MSSM id. 



Another scenario is given by the so-called cold electroweak baryogenesis where the elec- 
troweak phase transition is tachyonic and initiated at the end of (inverted) low-scale hybrid 
inflation. In this case, not the temperature dependence of the free energy but the inflation 



field turns the effective squared Higgs mass parameter negative 



12|-|l7l|. The electroweak 



phase transition occurs basically at zero temperature, and because of the spinodal instability 
induced on the Higgs field, its low momentum modes grow exponentially [la, ll8|, Q. This 
allows for a classical treatment of the dynamics. 

The generation of the baryon asymmetry in this scenario has been simulated on the 



lattice 



171], where the source of CP violation was assumed to be of the form 



HCP 

M2' 



tr {W^,W)^,) 



(1) 



where W denotes the SU{2)l field strength. A term of this form could originate in an 
effective action from a more fundamental theory at higher energies or from integrating 
out heavy fermions. In the present work we report on the simulation of cold electroweak 
baryogenesis using a CP-violating operator that was recently obtained by integrating out 
the quarks of the Standard Model |2o| and reads 



M2 



(2) 



Our main concern is if an operator of this form can successfully bias the Chern-Simons 
number and through the anomaly, baryon number, during tachyonic preheating and hence 
explain the observed baryon asymmetry via cold electroweak baryogenesis. 

It has often been argued in the literature that CP violation in the SM is too small to 
be able to generate the baryon asymmetry and we review this argument and the results of 
ref. 20| in the next section. In Sec. Illll we present numerical results of the lattice simulation 



of cold electroweak baryogenesis. Finally, in Sec. |IV]we conclude. 
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II. CP VIOLATION IN THE STANDARD MODEL 

It is often stated in the literature that the CP violation present in the SM is insufficient to 

explain the observed baryon asymmetry. These claims rest usually on the so-called Jarlskog 

I 

determinant 21| and we review this argument in the following. The basic observation is that 
physical observables cannot depend on the flavor basis chosen for the quarks; in particular 
transformations of the right-handed quarks leave the Lagrangian invariant since the weak 
interactions are chiral. Besides, the quark fields can be redefined absorbing one complex 
phase. The last fact implies that all CP-odd observables in the SM have to be proportional 
to 

J = SIS2S3C1C2C3 sm{6) = (3.0 ± 0.3) x 10"^ (3) 
with the Jarlskog invariant J given in terms of the Kobayashi-Maska wa p arametrization of 



the CKM matrix V with a CP- violating phase 6 as defined in refs. 2l|, |22|]. In addition, 
if two up- or down-type quark masses were degenerate, there would be no CP violation in 
the Standard Model since flavor basis transformation can in this case be used to remove the 
complex phase of the CKM matrix altogether from the Lagrangian. 

If one further assumes that the observable under consideration is polynomial in the quark 
masses, the simplest dimensionless expression that fulfills these constraints is found to be 
the Jarlskog determinant that has the form 

Acp = f ~^^Im Det |^m„mjj, m^mj^j 

= Jv-'' l[{m^, - ml^) llirhl - rhl^) ^ 10-^^ (4) 

where v is the Higgs vacuum expectation value and m^y^ denote the diagonalized mass 
matrices according to 

rridm^ = Dtti^D^ rrium^ = Um^JJ^ . (5) 

The identity in eq. (jl]) results then from the following relation of the CKM matrix (summa- 
tion over indices is only performed as explicitly shown) 



Im 



Va,VlV,,V}^ =JJ2 ^aceeuf, V = D . (6) 



According to this argument CP violation in the SM seems to be too small to explain the 
observed baryon asymmetry that is of order rj ~ lO^^*' and several proposals in the literature 
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aim at avoiding this bound. For example, it has been argued that at temperatures of the 
electroweak scale the CP violation might be only suppressed by the temperature rather 
than by the Higgs vev v as given in eq. (jlj), but this still is insufficient to be significant 
in a baryogenesis mechanism 23|, |2J| . In the context of electroweak baryogenesis coherent 



23 



2^ but finally 



scattering at the bubble wall has been suggested in the same work 
dismissed 25|. Furthermore construction of rephasing invariants containing derivatives of 



the Higgs field have proved to be incapable of being relevant in electroweak baryogenesis 
The approach in refs. 



20 



26|. 



27j is based on the effective action obtained by integrating 



out the quarks of the SM in the gradient expansion (more precisely an expansion in Lorentz 
indices, i.e. the covariant derivative expansion). The resulting effective action can potentially 
contain factors of the form 

These expressions vanish for degenerate quarks but are unlike the factors in the Jarlskog 



determinant not suppressed for small Yukawa couplings. The appearance of terms of this 



form was already argued in the context of cold electroweak baryogenesis in ref . |28| . Indeed 
when the quarks of the SM are integrated out, one finds in the effective action in next-to- 
leading order of the gradient expansion the operator [20 1 

2 / 

Z,W:^>W- {W+W- + WtW-) + c.c. ) , (8) 



J Kcp 



d X 



V 



8(47r)2 16 m2 
where J is the Jarlskog invariant given in eq. ([3]) and 

Kcp ~ 9.87. 



(9) 



As required, the operator is proportional to J and would vanish for degenerate quark masses. 
However, the latter information is hidden in the numerical coefficient kcp that is a function 
of the six quark masses ^ ^. The function is symmetric under the exchange of two families 



2C[ leads to no CP 



^ There is a recent claim 29] that an alternative method to the one used in ref. 
violation in the imaginary part of the Euclidean effective action in next-to-leading order. Moreover, using 
this method all operators considered by the authors (all CP-odd and several CP-even operators) of the 
imaginary part apparently vanish to this order in four dimensions. On the other hand, the next-to-leading 
order result in two dimensions (as presented in and confirmed with a different method in 30 1) does 



not vanish. A vanishing result in four dimensions seems implausible. We also stress that the results in [2o| 
are obtained on par with a large number of consistency checks and computed with the help of a computer 
algebra program^ 
^ The same work 



2^ also reported on CP violation in the real part of the Euclidean effective action that 
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but finite in the limit of one family becoming massless. The scale of the operator is hence 
given by the second heaviest family. 

This result was obtained in unitary gauge and the action can be rewritten in SU{2)l 
gauge invariant quantities. The charged gauge fields can be rewritten as 

^.^'H^ nr-^ilEE^l W-^^ W' = ^ (10) 



and similarly for the uncharged quantity 



Z, = -B,= ^ ^Ltl ' • (11) 



In the following we employ the operator (|8]) in a modified form in the lattice study of cold 
electroweak baryogenesis. This is necessary, since the gradient expansion used to obtain the 
result is only strictly valid in the case 

Z„ W^, « m, Z^,, W% « m\ (12) 

As mentioned earlier, the operator turns out to be finite for vanishing up and down quark 
masses and one finds that the next-to-leading order operator is suppressed by the charm 
mass as indicated in the notation of eq. (|H])- It is plausible that this is the appropriate mass 
scale that has to be used in the inequalities ( IT2|) . Still, the vev of Higgs field vanishes at 
certain points during the tachyonic phase transition, thus invalidating the expansion. 

In order to avoid that the dynamics is dominated by those points, we introduce a cutoff. 
The purpose of the cutoff is to suppress the CP violation in the region where the operator ^ 
is artificially large due to the break down of the gradient expansion. For the mechanism 
of cold electroweak baryogenesis, it is not essential that CP violation is active for Higgs 
configurations with almost vanishing vev. Two configurations with CP conjugate initial 
conditions will slowly drift apart, eventually leading to a difference in Higgs winding number, 
even if CP violation is only active in the regions of large Higgs vev. 

Instead of an expansion in gradients, an expansion in inverse gradients is feasable in the 
regions of small Higgs vev and higher order operators should come with a suppression of 
order 

A=«<^<%«^. (13) 



is however subdominant for cold electroweak baryogenesis due to parity conservation. 
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In practice, we insert A whenever the Higgs field appears in the denominator 

and c is just a constant fixed to satisfy the last condition. Even after introduction of the 
cutoff, the gradient expansion could be jeopardized by too large gauge fields. Means of 
justifying the gradient expansion could be to analyze the dependence of our results on the 
cutoff and the impact of the CP-violating operator a posteriori. 



III. NUMERICAL ANALYSIS 



Following the analysis of ref. |17l. I3ll-l34|. the SU(2)-Higgs model is discretized on a space- 
time lattice, and the classical equations of motion derived in a straightforward way. The 
implementation of the CP-violating operator hereby requires symmetrization in space and 
time, leading to implicit equations of motion in time, again in a similar fashion as for the 
operator ([1]) used in [l^]. The equations are solved by iteration and convergence of this 
iteration procedure imposes certain restrictions on the coefficient kcp, the timestep dt and 
the cutoff A, in order for the CP-violating force to not be too large. Ultimately, we are 
interested in the Higgs winding number and the Chern-Simons number that read 

= j dx\<P^d,<j)){<P^d,<P){<l)^dk<P), (15) 

Ncs = ^e,,fc j dx' (^WtW^, - \ea,cW^W'^wi^ . (16) 

where we used the temporal gauge, Wq = 0, and IV^^ denotes the corresponding field 
strength. 

The simulation generates random configurations ona.n^ = 64^ lattice, with lattice spacing 
arriH = 0.35. These are evolved using the equations of motion to mnt = 30, with timestep 
dt = 0.0125. Each configuration is run twice, using kqp = ±50, and subtracting the resulting 
Ncs and A^w values. This way, the ensemble of initial conditions is CP symmetric and the 
asymmetry is strictly zero at kcp = 0. This drastically reduces statistical noise, allowing 



a much clearer signal 33|. Correspondingly, only instances with integer difference in the 



Chern-Simons number between ±k,cp contribute. Statistics of this sort of observable is 
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FIG. 1: Two example configuration: The right configuration results in a difi^erence in the Higgs 
winding number for CP conjugate initial conditions, while the left does not. The Higgs field 
(straight curve) is normalized to v"^. The red and blue lines denote N^s and Nyj respectively. 

binomial such that the standard deviation is given by 



jumps 



22 7V, 



con figs 



2-^configs 



(17) 



\/-^configs T 

It will turn out that the asymmetry for all pairs of configurations is non-negative and that 
the result is overall not consistent with zero. 

As parameters we use mn = 2mw and a cut-off of cA^ = 50^ GeV^. The value of the 
cutoff corresponds to the estimate in ( fT3i) . In the simulation, we observe that the baryon 
asymmetry scales roughly with the fourth power of the cutoff what introduces a sizable 
uncertainty. A more sophisticated estimate of the cutoff and a more extensive presentation 
of the data will be given elsewhere (ssj . 

Figure [1] shows two example configurations, one of which contributes to the asymmetry. 
The Higgs field falls into the broken phase and performs a damped oscillation. In the 
meantime, energy is transferred into the gauge and Higgs fields, and their modes grow 
under influence of the CP violation. 

A first difference between the CP conjugate configurations occurs typically around trriH ~ 
6. At this time already a substantial amount of energy is transformed into the gauge fields 
and the Higgs vev is non-vanishing what is essential for the CP-violating operator ([2]) to 
be relevant. At the first minimum of the Higgs oscillation, where many zeros of the Higgs 
field are created j32| winding number is potentially generated. The presence of CP violation 
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FIG. 2: The final value of Higgs winding number for all 167 pairs of configurations. Black x-marks 
correspond to +kc'p, red +-marks to —kcp- In total 5 pairs of configurations contribute to a net 
baryon asymmetry. Deviations from integer values are rather small and due to lattice discretization 
errors. 

can hereby lead to a net baryon number between configurations with CP conjugate initial 
conditions. At late times winding number and Chern-Simons number agree and in the case 
at hand the Chern-Simons number follows the previously generated winding number. At 
the end, winding number approaches an integer valued vacuum while Chern-Simons number, 
containing thermal noise, oscillates for a longer time before it eventually settles into the same 
integer value |l7|. 

Notice that this behavior is quite different to the mechanism based on the operator ([T]). 
For a non-vanishing Higgs field this operator can be interpreted as a chemical potential for 
the Chern-Simons number. Accordingly, Chern-Simons number is already generated during 
the first roll-off of the Higgs field and in the first minimum of the Higgs field the winding 



33|. 



accommodates to the Chern-Simons number instead of vice-versa 

The final value of winding number is shown for an ensemble of 167 configurations in fig. [2l 
We note that all final values are very close to integers, showing that the lattice discretization 
errors are well under control. Only for one pair of configurations the winding number did 
not settle yet in an integer value at the end of the simulation. There are five configurations 
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FIG. 3: The average Higgs field (scaled by a factor 1/50), Chern-Simons number and winding 
number as functions of time. The average value 5/334 is superposed (dashed line) as well as the 
estimated error band (dotted lines). 



with a net winding number, indicated by arrows. Hence the asymmetry is 

5 



± 



(1/5 - 1/167) 



167 



0.015 X (1 ±0.44). 



(18) 



334 334 V 166 

Figure [3] shows the average of Higgs field, Chern-Simons number and Higgs winding 
numbers as functions of time. The average value 5/334 and the error estimate are indicated 
by blue dashed and dotted lines, respectively. The same picture arises as for the example 
configurations: Chern-Simons and winding numbers vanish until the first Higgs minimum, 
where winding number is generated potentially in a CP-violating way. The Chern-Simons 
number accommodates to the winding number subsequently. 

The baryon number density, reheating temperature and photon number density are given 

by 

ns = 3?2cs, Vo = — = ^^T^eh^ n-y = -—^T^^^. (19) 



4 30 ' 7 45 

For the Standard model g* = 86.25, and v = 246 GeV and we choose rriH = 2mw ~ 160 
GeV, what yields for the baryon asymmetry 



!^ = ^x(4.4±1.9)xl0-^ 



10 



(20) 



10 



which for the physical value kqp = 9.78 is four orders of magnitude larger than the observed 
asymmetry. 

IV. DISCUSSION AND CONCLUSIONS 

We presented first results of the simulation of cold electroweak baryogenesis utilizing the 
operator ([2]) that arises in the effective action of the SM from the CP violation in the CKM 
matrix. In order to make this operator applicable for small Higgs fields in the simulation, 
we introduced a cutoff. We chose a cutoff of cA^ = 50^ GeV^ and found a result that is 
four orders of magnitude larger than the observed baryon asymmetry. If the result would 
be dominated by the infrared modes of the Higgs field in the operator one would expect 
according to eqs. ffTOj) and f lT4|) a scaling as A~^^. Preliminary results of extensive computer 
simulations confirm that the result does rather scale as A~^ which we find very encouraging. 
A more detailed study of the dependence of the resulting asymmetry on the cutoff and the 



applicability of the operator (E]) in general will be topic of a subsequent work [35 1. 

We would like to comment on why Standard Model CP violation is operative in the 
proposed mechanism. Main requirement is that the quark masses enter in a non-polynomial 
way in order to avoid the Jarlskog determinant (jlj) as an upper bound on CP violation. In 
principle, this is easily achieved as can be seen in the Kaon system. However, baryogenesis 
typically takes place in a hot plasma and the temperature of the plasma provides a new 
energy scale. The temperature effects of the plasma render most processes in the infrared 



23 



24|. 



finite what makes it hard to avoid (j4|) or a similar bound depending on temperature 
Cold electroweak baryogenesis operates at zero temperature. 

In conclusion, we find it remarkable that a cosmologically viable baryon asymmetry could 
be created at a tachyonic electroweak transition, using only Standard Model CP violation. 



Although technical issues relating to the realization of low-scale inflation 
gradient expansion in this context persist 



and the 



20|], cold electroweak baryogenesis should be 



considered a serious candidate scenario for explaining the baryon asymmetry and deserves 
further investigation. 
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